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Abstract 

The A dependence of the production of hadrons in inelastic muon scattering and of 
the production of dimuons in high proton interactions are simply related. Feynman 
X distributions and z scaling distributions in nuclei are compared with energy loss 
models. Suggestions for new data analyses are presented. 
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I. Introduction 

When muons scatter from a proton or neutron the nucleon is excited and its struc- 
ture is changed because of the interaction of the muon with a quark in the proton. 
qq pairs which eventually become on-shell pions are produced and the yield of pions 
can be measured as a function of z, the ratio of the pion energy to the maximum pion 
energy allowable in the process. 

The time evolution of this state just after the interaction is unknown but it can be 
studied by measuring the same reaction in nuclei. In a sense the nucleus samples the 
emerging hadrons because the hadrons only rarely escape the nucleus without rescatter- 
ing. They impart energy to the spectator nucleons and hence lose energy themselves 
before coming on-shell. 

The final state (fully clothed) hadrons produced in this interaction are formed 
roughly at the dilated time (i//m)(r/i/c), where u is the energy imparted to the nucleon 
and r/i is the hadron size. For the data we study in this paper the energies, i/, are in 
the hundreds of GeV. 

Because the hadronization products cannot be calculated by perturbative QCD we 
carry out a study of models of the hadronization process in reactions involving muons. 

In one model the nucleon is excited and passes through the nucleus losing energy 
to the spectator nucleons in its path. This energy loss by the excited nucleon reduces 
the maximum energy available to the hadrons (mainly pions) it produces and will thus 
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affect the value of z = Eh/Emax- Because the muons are weakly interacting, they only 
interact at their full energy with a single target nucleon. (The evolution is shown in 
Fig. 1.) 

In a second model it is possible that the qq pairs, which eventually become the 
pions, are produced in the initial interaction. These pairs interact with the spectator 
nucleons on passing through the nucleus, imparting energy to them. (The evolution is 
shown in Fig. 2.) 

It is not difficult to estimate the effects of energy loss on either model. What 
is needed is an assumption about the interaction cross section. We can assume for 
simplicity that these off-shell hadrons have a cross section, ainei, for reinteraction that 
is the measured p-nucleon or pi-nucleon cross section, but we shall demonstrate how 
the z dependence depends on a range of ainei^s. We also need an estimate of the energy 
loss of the hadrons per collision. With this information the nuclear z distributions can 
be obtained. 

In a parallel study, also involving muons, but now in the final state, we examine 
the production of dimuons by protons: The incoming proton interacts with spectator 
nucleons and makes soft minimum bias collisions before making the high Q"^ dimuon 
pair. Thus the incoming proton's energy can be decreased, thereby reducing the energy 
available for producing the rare dimuon pair. This shifts events to lower Feynman x. 
X = E(dimuon)/Einax(diinuon). (The evolution for this reaction is shown in Fig. 3.) 

In muon production of hadrons the A dependence is often revealed by exam- 
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ining R(z) = Na(z)/Nd(z) while in dimuon production by protons the quantity is 
R(x) = Na(x)/Nd(x). N is the number of events at either x or z. The comparison 
is usually made with the deuteron, D, to include possible p-n differences. 

A most interesting value of R is the limiting value at z or x equal to unity. In the 
hadron production case, if the the muon interacts with the last nucleon in its path, 
there will be no energy loss. Only such events will appear at z=l. In the dimuon case, 
the events at x=l arise from those interactions where the incident proton made the 
dimuon on the first nucleon in its path, the weakly interacting dimuon pair passing 
unscathed through the nucleus. If a„ is the (Glauber) probability of making n collisions 
in a proton-nucleon minimum bias interaction in a nucleus, nUn is the probability of 
making a rare high collision like a muon scatter or the production of a dimuon 
pair or vector meson such as a J/^ . Since the probability of making one collision 
is then na„/n = a„ the fraction of times the collision takes place in a first or last 
collision is just X^a^/ Z)na„ = V < n > Thus these endpoints depend only on Glauber 
probabilities and not on the detailed mechanisms. This general argument shows that 
R cannot be a constant, independent of z or x. It is not difficult to calculate these 
endpoint values using a Woods-Saxon spatial distribution for the nucleons and the 
total inelastic cross section for the nucleon-nucleon scattering. That endpoint does not 
depend on the energy loss function or its magnitude. 

Actually, we really need to know the inelastic cross section for an off-shell hadron 
on a ground state nucleon if we wish to calculate R at any other value than R(l). 
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II. Hadron Production in Inelastic Muon Interactions 

Figure 4 shows a plot of R(z) taken from the data of ref. [|T| The authors have 
chosen to separate their data sample into a low - low xbj bin and a high - high 
XBj bin, where is the four momentum transfer to the nucleon and XBj is Bjorken 
X. As these authors have shown, and because it can be seen that the measured points 
are almost identical in the two samples, there is apparently no nuclear dependence in 
the data on these variables. As a result, and to improve the statistical accuracy of 
our comparison, we have suitably averaged the two results which are shown plotted in 
Figure 5. 

To examine the relative merits of fits of the data to the conclusion in ref. |]l|that 
there is no A dependence, we have examined the relative fo^^ ^ linear fit hypothesis 
to the data shown in Figure 5. For the fit to R = 1, the for the five measured points 
is 7.2. It is 3.0 for the best linear energy loss fit shown in the figure. (Eliminating the 
low statistics high z point hardly changes the slope of the falloff but improves the 
by 2 rather than 1 unit.) The authors of ref. remark that there are large systematic 
errors in the lowest z point. Without that point we find that the x? is .53 for the linear 
fit and 3.2 for a flat fit with no A dependence, i.e. R = 1. Thus an A dependent effect 
seems to provide a better fit to the data, suggesting that there is a nuclear effect to be 
understood. 

In several previous papers we have examined the effects of energy loss in nuclear 
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interactions. 0] The calculation of R(z) is simple if one knows the form of the 
energy loss per collision. Various models ^ have been proposed recently, giving 
different assumptions for the a/s dependence of the energy loss. We make use of the 
model we used in 1987 since it was demonstrated to be in agreement with the 
data observed in low pt production and was the functional form for the energy loss 
per nucleon obtained by examining the energy loss found in the ISAJET minimum 
bias model of F. Paige[|^, widely used by particle physicists. That energy loss is 
approximately given by d^jdn = constant = (3, where n is the number of collisions. 
Thus in our model the energy loss varies as the square root of the laboratory energy, 
\/E This energy dependence lies between the values and of refs. and 
respectively. 

We first consider the model of the struck proton hadronizing outside the nucleus 
into pions. In Fig. 6 we have plotted our results for several values of z/ = — E^/. 
After these calculations were made we obtained the actual u distributions obtained in 
the experiment. That average value is z/ = 170 GeV which one can see from Figure 
6 would fall well on the experimental results. All the theoretical curves predict a rise 
above R = 1 at low z due to the sliding back of events to lower z due to energy loss. 
While this appears in the data, the authors of ref. 1 caution that there are large 
systematic errors in the lowest z point. Note that the point at z = .68 falls well below 
R = 1. For this calculation we used d^/s/dn = .2 GeV, which we had found in earlier 
work fit the available data on dimuon and J/\E' production [0. The theoretical curves 
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have been corrected to take into account the bin size used in the data. The N(z) data 
for both Xe and D were fitted to a sum of two exponentials of the form e~"^ to get the 
best representation of the input to the energy loss calculation. To show the sensitivity 
of the calculations to both the inelastic cross section, aimi, and d^/s/dn , we show in 
Figures 7 and 8 how variation in these parameters affects the results. These plots also 
show the asymptotic limits. There is roughly a trade-off of d^/s/dn of .1 GeV for a 
change in the inelastic cross section of 10 mbs. 

We now turn to the second calculation, namely the model in which the qq pairs are 
formed at the collision and lose energy on the way out of the nucleus. Fig. 9 shows our 
results for an inelastic pion-nucleon cross section of 20 mb and an energy loss d^/s/dn 
= .2 for different values of the energy of the pion. The dependence on energy is small 
and all curves show a depression below R = 1. There is a good fit to the data if we 
again omit the lowest z point as suggested in ref. Figure 10 shows the effect of 
varying the energy loss when u = 170 GeV and atot = 20 mbs. 

Fig 11 shows the effect of varying the cross section of the pion using the average 
value of u of 170 GeV and the same energy loss parameter. Here we see that 20 mb 
appears to give the better fit than 10 or 31 mb. Once again we note different limits at 
R = 1. 

Comparing Figs 9, 10, and 11 with Figs 5, 6 and 7, we note that in the case of the 
pion energy loss the slopes of R vs z are smaller and the main effect is a depression of 
R below unity. 
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We conclude that while the data are not very precise they do not rule out the 
presence of energy loss mechanisms, as suggested by W. Buszaf^. 

III. Dimuon Production in p-A Collisions 

Dimuon production in p-A collisions can be calculated on the same energy loss 
model, but is slightly more complicated since the dimuon production cross section is 
energy dependent and the shape of the cross section depends on the invariant mass of 
the dimuon pair. 

The effect on R{xf ) = <Jh-A(yXF)/Aah-p{xF) of any loss in energy comes about 
because of the dependence of the Drell-Yan cross section which varies as e^^^l^. 
M is the invariant mass of the muon pair. The energy loss therefore produces two 
effects a) a reduction in the dimuon yield and b) a displacement of events to lower xp. 
Before demonstrating the effects of this energy loss on R{xp ) we present a "back of 
the envelope" calculation of R = , which shows the qualitative features. 

l)i?=l-(7)(M/v/i)(/3/v/i)((n)-l) 

This formula demonstrates how the various parameters enter into the dimuon "de- 
pression". ( Empirically the last bracket in eq. 1) varies very roughly as In A so, for 
small departures from unity, R will vary as A raised to a small constant.) While R 
should approach unity at large s, it will not be unity at lab energies as high as 800 
GeV, nor will R{xf ) be independent oi xp . 
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We now turn to two pieces of data that illustrate the effect of energy loss. We 
use the Drell-Yan formalism and the Duke-Owens form factors for the calculation ||1CI||, 



rather than using the empirical values for the energy dependence of the cross section, 
since we have checked that they give essentially the same results. 

Fig. 12 shows the data of the NAIO collaboration on dimuon production by 140 
and 286 GeV pions[|n|]. Superimposed on their data are our calculations for .2 and .4 
GeV energy loss. There appears to be a clear M dependence as well as deviations from 
R=l in the xp distributions. The data is not very precise but the general features of 
energy loss are borne out. 



Fig. 13 shows the 800 GeV data|T^ of Fermilab expt. E 772. We have obtained 



the unpublished mass distributions |^ for each xp bin so we can make a comparison 
between the data for different values of dy/s/dn . The published xp distributions for 
their W to Deuterium ratios show a deviation in R from unity at large xp . Our 
calculations for d^jdn = 0.2 are superimposed on the data. 

We conclude that there are A dependent effects in dimuon production which can 
be accounted for using the same energy loss mechanism that can account for the A 
dependence in inelastic muon scattering and using the same approximate energy loss 
parameter. 

IV. Discussion 

Unfortunately, none of the data extend to very large values of z or x to enable the 
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asymptotic values at R(l) to be compared with the Glauber prediction, verifying the 
importance of the energy loss in the most unambiguous way. There are other recent 
estimates for the functional form of the energy loss P but the present data cannot 
easily be used to discriminate among the various models. 

However the question of the form of the energy dependence of d^fsjdn can be more 
easily examined with present data on hadron production by muons. E665 can use all its 
data not just their published low and high - Bjorken x bins to improve the statistical 
accuracy for a new study: the separation of the present E 665 muon inelastic data into 
separate v bins to study R(z/). Newer E665 data could be analyzed to study R for 
a wide range of A. For example, with our parameterization dy/s/dn , the laboratory 
energy loss, dE/dn, would vary as \fE so that the loss parameter would clearly change 
in the region covered by the E665 data (110 to 490 GeV). While it is difficult to study 
the energy dependence of d\/s/dn in dimuon production, since different accelerator 
energies are needed, it is simple to select v from the inelastic muon data to make the 
same study. 

This work extends to muon inelastic scattering and dimuon production the need 
for energy loss in traversing nuclear matter in accounting for A dependences which we 
showed earlier could account for various A dependences in minimum bias reactions. 
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Figures 

Fig. 1 Hadron Production by Muons: An incident muon interacts with a nucleon 
which rescatters off another nucleon in the nucleus. The nucleon struck by the muon 
hadronizes outside the nucleus. 

Fig. 2 Pion Production by Muons: An incident muon interacts with a nucleon 
producing an off-shell pion which rescatters off nucleons in the nucleus 

Fig. 3 Dimuon Production by Hadrons: An incident proton makes soft non- 
perturbative scatters before making a high collision with another nucleon producing 
an excited proton which hadronizes outside the nucleus. 

Fig. 4 The z dependence of the ratio of hadron production in Xenon and Deuterium, 
R(Xe/D), for high (open circles) and low (solid circles) and xsj- The data are taken 
from ref. 1. Table XXXI. 

Fig. 5 R vs z for the combined sample of high and low data. The solid curve is 
a least-square fit to a straight hue fit to the data. 

Fig. 6 Hadron Energy Loss: R vs z for the combined sample of high and low 
data. Theoretical curves are shown for three values of the energy transfer, v. The 
mean value of v for the data is 170 GeV. 

Fig. 7 Hadron Energy Loss: 7 Theoretical calculations of R for different values of 
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d^/dn . The inelastic p-p cross section is set at 31 mb. The value of R at z=l is set 
by this cross section and the Woods-Saxon nucleon spatial distribution. 

Fig. 8 Hadron Energy Loss: Calculations of R for different inelastic cross sections. 
The value of d^/dn = .2 GeV. Note the different asymptotes. 

Fig. 9 "Pion Energy Loss" : R vs z for the combined sample of high and low 
data. Theoretical curves are shown for three values of the energy transfer, u. The 
pion-nucleon cross section is taken as 20 mb. d^/s/dn = .2. The mean value of u for 
the data is 170 GeV. 

Fig. 10 "Pion" Energy Loss: Theoretical calculations of R for different values of 
d^/s/dn = .1, .2, .3, and .4. The inelastic pi-p cross section is set at 20 mb. The 
value of R at z=l is set by this cross section and the Woods-Saxon nucleon spatial 
distribution. 

Fig. 11 "Pion" Energy Loss: Calculations of R for different inelastic cross sections, 
10, 20 and 31 mb. The value of dy/s/dn = .2 GeV. u = 170 GeV. Note the different 
asymptotes. 



Fig. 12 Data of ref. ||Tl| for dimuon production by pions. The predicted falloff of R 
with the dimuon invariant mass, M, as well as with Feynman X are shown along with 
the theoretical energy loss predictions. 



Fig. 13 Dimuon production by 800 GeV protons. R vs a;^ . Data from ref.[|T2 
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